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Abstract. We consider Glauber dynamics of classical spin systems of Ising type 
in the limit when the temperature tends to zero in finite volume. We show that 
information on the structure of the most profound minima and the connecting 

i-rt I saddle points of the Hamiltonian can be translated into sharp estimates on the 

P3 ■ distribution of the times of metastable transitions between such minima as well 

O I as the low lying spectrum of the generator. In contrast with earlier results on 

such problems, where only the asymptotics of the exponential rates is obtained, 
we compute the precise pre-factors up to multiplicative errors that tend to 1 as 
T I 0. As an example we treat the nearest neighbor Ising model on the 2 and 
3 dimensional square lattice. Our results improve considerably earlier estimates 
obtained by Neves- Schonmann [NS] and Ben Arous-Cerf [BC] and Alonso-Cerf 

Xf-\ I [AC]. Our results employ the methods introduced by Bovier, Eckhoff, Gayrard, 

t^ ! and Klein in [BEGK1,BEGK2]. 
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1. Introduction. 

c^ I Controlling the transitions from metastable states to equilibrium in the stochastic 

S ' dynamics of lattice spin systems at low temperatures has been and still is a subject 

i-rt ■ of considerable interest in statistical mechanics. The first mathematically rigorous 

^ ! results can be traced back to the work of Cassandro et al. [CGOV] that initiated 

the so-called "path-wise approach" to metastability. For a good review of the earher 
literature, see in particular [Va]. All the mathematical investigations in the subject 
require some 'small parameter' that effectively makes the timescales for metastable 
rS I phenomena 'large'. The somewhat simplest of these limiting situations is the case 

c^ I when a system in a finite volume A C Z*^ is studied for small values of the tempera- 

ture T = 1/(3. In systems with discrete spin space one is then in the situation where 
the dynamics can be considered as a small perturbation of a deterministic process, 
a situation very similar to what Freidlin and Wentzell [FW] called 'Markov chains 
with exponentially small transition probabilities'. Consequently, most of the work 
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2 A. BOVIER AND F. MANZO 

concerning this situation [0S1,0S2,CC,BC,N,NS] can be seen as extensions and im- 
provements of the large deviation approach initiated by Freidhn and WentzelL This 
consists essentially in identifying the most likely path (in the sense of a sequence 
of transitions) and proving a large deviation principle on path-space. While this 
approach establishes very detailed information on e.g. the typical exit paths from 
metastable states, the use of large deviations methods entails a rather limited preci- 
sion. Results for e.g. exit times r are therefore typically of the following type: For 
any e > 0, 

P (e/3(A-.) < ^ < ^/3(A+.)) t 1, as /5 r oo 

where A can be computed explicitly. Similarly, one has results on eigenvalues of the 
generator that are of the form 

lim/5-MnA,(/3) = 7^ 

/3Too 

with explicit expressions for the 7^ (see e.g. [FW,S]). From many points of view, the 
precision of such results is not satisfactory, and rather than just exponential rates, one 
would in many situations like to have precise expressions that also provide the precise 
pre-factors. This is particularly important if one wants to understand the dynamics 
of systems with a very complex structure of metastable states, and in particular 
disordered systems. (For a rather dramatic illustration, see e.g. [BBG1,BBG2] where 
aging phenomena in the random energy model are studied.) Another drawback of 
the large deviation methods employed is that they are rather heavy handed and 
require a very detailed knowledge of the entire energy landscape, a requirement that 
frequently cannot be met. 

In two recent papers [BEGK1,BEGK2] a somewhat new approach to the problem 
of metastability has been initiated aiming at improving the precision of the results 
while reducing at the same time the amount of information necessary to analyze a 
given model. To achieve this goal, the attempt to construct the precise exit paths is 
largely abandoned, as are, to a very large extent, large deviation methods. 
The general structure of this approach is as follows. In [BEGK2] the notion of a set 
of metastable points is introduced. The definition of this set employs only one type 
of objects, namely Newtonian capacities (which may also be interpreted as escape 
probabilities) . If such a set of metastable points can be identified, [BEGK2] provides 
a general theorem that yields precise asymptotic formula for the mean exit time 
from each metastable state, shows that this time is asymptotically exponentially 
distributed (in a strong sense), and states that each mean exit time is the inverse of 
one small eigenvalue of the generator. While [BEGK2] assumes reversibility of the 
dynamics, in [E] is shown that almost the same results can be obtained in the general 
case. Thus, the analysis of metastability is essentially reduced to the computation 
of Newtonian capacities. The great advantage of such a result is that capacities 
are particularly easy to estimate, due to the fact that they verify a particularly 
manageable variational principle. This fact is well known and has been exploited 
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in the analysis of transience versus recurrence properties of Markov chains (see e.g. 
[DS]); however, its particular usefulness in the context of metastability seems to have 
been noticed only in [BEGKl] where it was used in the context of reversible discrete 
diffusion processes motivated from certain mean field spin systems. 
In this short paper we will show that the approach is even more efficient and simple 
in the context of the zero temperature limit of Glauber dynamics of spin systems in 
finite volume. We will show that in rather general situations, capacities in this limit 
can be computed virtually exactly in terms of properties of the energy landscape, and 
therefore all interesting properties of the dynamics can be inferred from a (not overly 
detailed) analysis of the energy landscape generated by the Hamiltonian considered. 
As a particular application that should illustrate the power of our approach, we apply 
the general results to the Ising model (in two and three dimension). 

2. The general setting and the main theorem. 

In this section we set up the general context to which our results will apply. It will 
be obvious that Glauber dynamics of finite volume spin systems at low temperatures 
provide particular examples. We will consider Markov processes on a finite state 
space VL (the configuration space). To define the dynamics, we need the following 
further objects. 

1. A connected graph Q on VL. We denote by E{Q) the set of the edges in Q. 

2. A Hamiltonian H -.Vt -^M. also called energy. 

3. The Gibbs measure Q{x) := ^exp{—f3H{x)), where Z is the normalization 
factor called partition function, and (3 is the inverse temperature. 

We consider transition probabilities P{x,y) such that if {x,y} G E{Q), P{x,y) > 0, 
and P{x, y) = ii X y^ y and {x, y} ^ E{Q). We assume moreover that the transition 
probabilities are reversible with respect to the Gibbs measure, i.e. 

Q{x)P{x,y)=Q{y)P{y,x). (2.1) 

We will also make the simplifying assumption that any existing transition in the 
graph is reasonably strong, i.e. we assume that there exists a constant C > such 
thatfl 

P(x, y) + P(y, x) > C V {x, y} e E{g), (2.2) 



by reversibility, (p.2|) is equivalent to 

Pix, y) > — , ,,^, . ^, .^. V {x, y} e E{g). (2.3) 

1 + exp(-/3(i7(|/) - H{x))) 

To be able to state our results we need some further notations. 



^The constant C will typically be of the order of the inverse of the maximum coordination number 
of the graph Q. 
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1. Given a one- dimensional subgraph cj, we write u; : x — *> / if the subgraph 
has one end in x and the other end in /. One dimensional subgraphs have a 
natural parameterization ujq, . . . , ujk, where K := |cj| — 1, V/c = 0, . . . , fi' — 1 
q{ujk,^k+i) > and u : ujo ^ ujr- 

2. Let H{{uj}) := raaxz^uj H{z). For x & Q and / C fi, we introduce the commu- 
nication height, H{x,I), between x and / as 

H{x,I) := min H{{lu}). (2.4) 

Moreover we define the set of saddle points for x and / by 

<Sx,i '■= \ z E Vl ; 3 CO : X ^ I with z E cu and H{z) = H{x, I) > (2-5) 

3. Furthermore, we define the set of points 

Di := {z ; H{S,,^) < H{S,j)} (2.6) 

These will be the points that are 'closer' to x than to /. 

4. For any set A C fi, we define its outer boundary dA as the set of all points in 
A from which an edge of Q leads to its complement, A'^. 

5. For z E dDl, let p^ := Y^y'eDi ^i^^v') and p^ := Y.x'en\Di P{z,x'). We let 

Cx,i ■■= V J4^. (2.7) 

6. Let W^ := {y ; H{y) < H{x)}. For x E Q, we set r(x) := H{S^^wJ - H{x). 
If X is not a local minimum of the Hamiltonian, r(a;) = 0. If x is a global 
minimum of the Hamiltonian, we set T{x) = oo. 

7. For the process Xf starting at x, we define the hitting time to the set / C ^2 as 
Tf := inf{0 <t^N; XtEl}. _ 

8. We denote by M. the set of all local minima of H. We call a subset A1 C A^ a 
set of metastatic states, if a point that realizes the absolute minimum of H is 
contained in A4 and if, for all y E A^\A1, T{y) < mina;g_A4 T{x). It is important 
to realize that for given H, Ai may often be chosen in different ways. The idea 
will be that we will observe the process only at its visits to Ai. Thus, the actual 
choice of Ai will depend on how much information we want to retain about the 
detailed behavior of the process. Note that this definition implies that for all 
z^M, H{z,M) >H{z,W,). 

9. Finally, for x E Ai, we need to define the quantity 

N^.-^fzEn-. {H{z) = H{x)}n{H{x,z) <H{x,M\x)}\, 

which represents the degeneracy of the minima of the Hamiltonian. 

We can now formulate the main general result of this paper in the general setting. 
Let us consider some set At of metastable points. To be able to formulate concise 
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and general results, we make some further assumptions that will be true for "generic" 
Hamiltonians.0 

(hi) For any x i^y E M, T{x) ^ T{y). 

(h2) For any x ^ y E M., S^^y consists of isolated single points^. 

Theorem 2.1. Let M. he a set of metastable states for the Hamiltonian H satisfying 
the conditions (hi) and (h2) above. Forx e M., set M.^ := {y g ^Vf ; H{y) < H{x)} = 
M. n Wx- Let r(x) := t^^, S := Sx,Wx ^'^^ ^x '■= C^yv^- Then there exists 5 > 0, 
independent of (3, such that for any x E M., 

(^) 

E t{x) = Ar,C-ie^^(")(l + o{e-^^)) (2.8) 

(ii) there exists an eigenvalue X^ of 1 — P such that 

K = ^^Al + o{e-^')) (2.9) 

E t{xj 

(Hi) if (px is the right- eigenvector of P corresponding to Xx, normalized so that 
(f)x{x) = 1, then 

0.(y) = p(r,^<r|^J+o(e-^^)) (2.10) 

(iv) 

P (r(x) > t E r(x)) = e-'^^+°^''''^\l + 0(6"'^^)) (2.11) 

In Section ^ we will apply Theorem 2.1 to a well known situation, the kinetic Ising 
model, in the limit of vanishing temperature. 

Let us anticipate our main result about the kinetic Ising model, referring to Section 
^ for precise definitions and notation. 

Theorem 2.2. Consider the kinetic Ising model with Metropolis dynamics in dimen- 
sion d = 2ord = 3ina torus A'^{1) with diameter I. The magnetic field < h < 1 
is chosen such that 2{d — l)/h is not an integer. Then, 

the two configurations —1 (all minus spins) and +1 (all plus spins) form a metastable 
set (if the magnetic field h is positive, +1 is stable), and 

• In dimension 2, let £2 '■= [f 1 o,nd T2 '■= 4^2 — /i (-^2 ~ ^2 + 1) be the diameter 
and the activation energy of the "critical droplet", respectively. Then, 

Er(-l) = l^—ef'^^l + o(e-^^)) = ^hef'^'il + oik) + o{e-^')) (2.12) 

O to — 1 lo 



^Note that for any Hamiltonian, one may select different sets of metastable points. The require- 
ments (hi) and (h2) depend on the Hamiltonian as well as on the choice of A4. 

^ In the appendix we will explain how one can proceed to obtain comparable results in the case 
when this condition is not satisfied. 
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• In dimension 3, let £3 := [|] and 

Ts := {6il - 4^3 + 4^2) -h{il-£l + el-£2 + l) 

be the diameter and the activation energy of the "critical droplet", respectively. 
Then, 

(2.13) 

iifere as in Theorem 2.1, 5 > Q is independent of 13 (hut depends on arithmetic 
properties of h). 

Remark. Note that in our model we flip at most one spin per time step. In 
continuous time dynamics the mean transition times would be lowered by a factor 

1/IAI. 



The above Theorem shows how the results of Theorem 2J. can be applied (via the 
analysis of the energy landscape carried out for the Ising model in [NS] and [AC,BC]) 
to the so-called Freidlin-Wentzell regime. Notice that the methods of [BEGK2] can 
be applied in a very similar way to situations where the volume grows with (3 to com- 
pute "exactly" the probability of first appearance of a critical droplet (a preliminary 
problem for the infinite- volume metastability carried out in [DeSc]). 

3. Basic tools. 



Theorem 2J. relies on Theorem 1.3 in [BEGK2] that links relative capacities of 

metastable sets to mean exit times and to the low lying spectrum of 1 — P. The 

additional work needed to prove Theorem |2.1| will be to estimate capacities in terms 

of the Hamiltonian H, and to show that the hypotheses of Theorem 1.3 in [BEGK2] 

are satisfied in our setting. 

Let us state Theorem 1.3 in [BEGK2] specialized to our case. 

In their context, a set A^ G f2 is called a set of metastable points in the sense of 

[BEGK2] if 

The set M. is generic in the sense of [BEGK2] if for any x,y E A4, and / C A4, 
-^ — — tends either to zero or to infinity, as /5 t 00, and if the absolute minimum 



of the Hamiltonian is not degenerate. 



Theorem 3.1. (Theorem 1.3 in [BEGK2]) Let M. he a generic set of metastahle 
states in the sense of [BEGK2], and let for x e M, M.x o-nd t{x) he defined as in 



Theorem \2. 4 Then, for any x E M., the following holds: 
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(i) 

N 

(a) for any x E M., there exists an eigenvalue X^ of 1 — P such that 

moreover, the eigenvalues of 1 — P not corresponding to any x E Ai are in the 
interval {c\Q\~^ inf^gn P (r_^ < r^) , 2] for some positive constant c. 
(Hi) if (px is the right- eigenvector of P corresponding to X^, normalized so that 
(f)x{x) = 1, then 

My) = ^{r!<rl,J+oil) (3.4) 

(iv) for any x E Ai, for any t > 0, 

P {t{x) > t E t{x)) = e-*(^+°«) (1 + o(l)). (3.5) 

Here o(l) stands for a small error that depends only on the small parameters intro- 
duced via ( p.l|) and the non-degeneracy condition following it. 

We leave it to the reader to verify that this theorem is indeed a special case of the 

more general result stated in [BEGK2]. 

Theorem 12.11 will follow from Theorem 13.11 since in the finite-volume and (3 —>■ oo 



regime, we compute P (tj^^ < t^) and show that local minima of the Hamiltonian 
are metastable states giving at the same time the value of the nucleation rate in the 
limit /? ^ oo. 
The key estimate is the following Lemma. 

Lemma 3.2. \/x,y G M. such that S^^y is a set of isolated single points, 

P [t^ < r^) = C,^j,e-^W'5--)-^("»(l + o{e-^^)). (3.6) 

We will explain in the appendix how our method can be extended to situations where 
the saddles are degenerate. In this case the pre-factor C^^y does not have the nice 
form in (p.7|) but can still be computed explicitly in terms of small " local variational 
problems" . 

Lemma 3.3. Let x he a minimum for the Hamiltonian. Then, x is a metastable 
state (in the sense of [BEGK2]) in the set Ai := {y ; T{y) > r{x)}. 

Clearly, Theorem |2.1| immediately follows from Theorem |3.1|, Lemma p.2| and Lemma 
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4. Proof of Lemmata ^T2] and ^ 



In order to prove Lemmata |3]^ and 

[BEGKl]. 

The following Lemma corresponds to Theorem 6.1 in [Li]. 

Lemma 4.1. (Dirichlet representation). 

Let n^:={h:n^ [0, 1] ; h{x) = 0, h{y) = 1} and 



we make use of many ideas contained in 



//\l2 



h{x")] 



x',x"en 



Then, 



-l3H(x) 1 



z 



(4.1) 



(4.2) 



n 



Proof. See [Li], Chapter II. 6. 

Note that the left-hand side of ( [4. 21) has the potential-theoretic interpretation of the 
Newtonian capacity of the point y relative to x (i.e. the electric charge induced on 
the grounded site x when the potential is set to 1 on the site y) . The Dirichlet form 
is just the electric energy, and the minimizer h* is the equilibrium potential, with the 
probabihstic interpretation h*{z) =F (vy < r^). 

The strength of this variational representation comes from the monotonicity of the 
Dirichlet form in the variables P{x',x"), expressed in the next Lemma, known as 
Rayleigh's short-cut rule (see Lemma 2.2 in [BEGKl]): 

Lemma 4.2. Let A be a subgraph ofQ and let Pa denote the law of the Markov chain 
with transition rates, for u ^ v, defined by Pa{u,v) := P{u,v)I{{u,v} G E{A)}. If 
X and y are vertices in A, then 



P {t^ < r:) > Pa (r; < t:) 



(4.3) 



Proof. The proof follows directly from Lemma |4.1| and can be found in [BEGKl] . D 



The following Lemma corresponds to Lemma 2.3 in [BEGKl] and is a well known 
fact (see e.g. [DS]). 



Lemma 4.3. (The one dimensional case). 

Let oj be a one- dimensional subgraph ofQ, K := \uj\ — 1 and let {ujn}n '■ {0, 

Q be such that Vn < K, q{ujn,ujn-i) > 

-1 



K} 



P (t'^° < t'^") 



'^ p-li(H{uQ)-H{uj^)) 
.n=0 



F'y^m^n+i) 



(4.4) 
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Remark: Lemmata |4.2| , |4.3| and ( p.2|) immediately give the following bound: Vx', / 
s.t. 5a;/,/ is made of simple points, 

-1 



p (rf < t:: 



> 



■^-1 g-/3(H(x')-H{c^„)) 



E 



n=0 



P {uJn,UJn+lj 



> 



> c 



,n=0 



g-/3(H(a;')-//(a;„)) _^ g-/3(/f(x')-/^(a;„+i)) 



> £g-/3(/^(5,/,,)-H(x')) ^^ _ g-/35\ 



(4.5) 



for any choice of the subgraph u : x' ^ I having its maximum energy in S^'j- The 
constant C is the same as in 



Proof of Lemma ^7^. Let F := H{Sx,y) — H{x). 
We consider the surface Z := dD^. Notice that: 

1. i_) :^ ^x y ^ ^ 

2. 3<5 > such that \fzeZ\S, H{z) > H{S) + 5. 

3. Z is the outer boundary of a connected set that contains x. 

Remark: In what follows, any other surface with properties 1, 2, and 3 would give 
the bounds we need for the proof. The quantity C'^y defined with respect to the new 
surface differs from Cx^y by a factor 1 + o{e~^^). 
We set Dx := D^, Dy':=Q\{ZU Dx), Z- := dZ n Dx and Z+ := dZ n Dy. 

(1) The upper bound. 

We use Lemma W^ with h{x') := if x' G Dx and h{y') := 1 ii y' E Dy] we choose 
h{z) for z E Z in an optimal way. In the rest of the space we choose h{x) = 1. 
We have 

P {r^ < r:) < ^^Hh) = 
= J2 e-^(^(^)-^("))p(2, x) {p,h\z) + p.(l - h{z)f) + o(e-^^), (4.6) 

z&Z 

where we used reversibility. The small error comes from the mismatch on the bound- 
ary of Dx that lies higher than the saddle hight. 

The quadratic form ph"^ + p(l — hY has a minimum for h = ■^. Hence, we can 
saturate the inequality ( [4.6|) and get 

(Lh.s. of (pD) < C,,,e-^^ (1 + e~^') . (4.7) 

(2) The lower bound. 

We consider the subgraph A obtained by cutting all the connections to the vertices 

in Z\S. 

We use Lemma E^ to bound the original process by the restricted process. 



10 



A. BOVIER AND F. MANZO 



We use (|4.5|) to estimate the probability to reach x' G Z^ and the probabihty to go 

from y' G Z^ to y: 

By the strong Markov property at time rf, we have 



PAh^'<r:) = 5^ Pa (r; < r|, J Pa (r;<r:) 

= e-^r ^ Pa (r,^<r|) Pa (r^^<r|) 



(4.8) 



zes 



where we used reversibihty. 

Now, 



Pa (r: < rj) = J] P {z, x') Pa (rf < r^') . 



(4.9) 



x'&Z- 



We bound the last factor using a standard renewal argument (see e.g. [BEGKl] 
Corollary 1.6) that yields if z' G D^ the last term is exponentially close to 1: 



Pa r| <rr 



A \Ts 



< 



Pa {r^is < <') ~ 
\S\e-^^ 



Pa (rr < r^;) 



< e 



-/35 



Ce-^(r-5')(i_e-/35'^ 
where we used (^75|). By putting together (^]9|) and ( [4.10|) we get 

Pa [rl < rJ) > p, (l - e-^-^) 
We use the same procedure to bound the last term in ( |4.8| ): 

Pa (r; < r:) > $^ P (^, y') Pa (r,^' < rf ) . 

y'&Z+ 



Again, the same arguments leading to ( [4.10|) show that the last term in this 
exponentially close to 1: 

>^<-f3{H{z)-H(y')) 



^A \^l^ ^ 'y ) - ^-p(H(Sy^^,)-H(y')) ^^ _ ^^pS' 

We put together ( |4.12|) and ( [4.13| ) and get 



< e 



-ps 



Pa (r; < r^) > 



1-e- 



Going back to (^, we get from (|4Tl| ), (gTT^) 

Pa (r,^' < r:) > C.,,e-^^ (l - e'^^ 



(4.10) 

(4.11) 

(4.12) 
sum is 

(4.13) 

(4.14) 

(4.15) 
D 
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Proof of Lemma [g. Sj . For any z ^ Ai v/e know that by definition of A^, we liave 
tliat T{z) < min^^gTn r(a;) = F. In view of Lemma |4.1| and tlie lower bound (4.5) 
we only need to show that this implies that H{z,M) — H{z) < F. Now let u ^ Ai 
be the point that realizes the minimum of the energy among the states such that 
H{z,u) < H{z,M.). For such a point, by definition, H{u,M.) — H{u) = F(-u) < F. 
But clearly, H{z,M.) — H{z) < H{u, Wu) — H{u) < F(x), and we are done. D 

5. The Ising case. 



In this section we want to illustrate the strength of Theorem p.l| in a well known 
context, namely the stochastic Ising model on the (i-dimensional lattice. In this case 
the state space is ^2 = { — 1, +1}^, where A = A(L) is a torus in 7/ with side-length 
L. For a &VL, the Hamiltonian is then given by 

H{a) = H^{a) = -]^ J2 a,a,-/ij]a„ (5.1) 

<i,j>€A ieA 

where the first sum concerns all the pairs of nearest neighbor sites in A. 
Let a'^ be the configuration that differs from a only in the value of the spin of site i and 
[a]+ denote the positive part of the real number a. We will consider for definiteness 
only the case of the Metropolis dynamics, i.e. the transition probabilities are chosen 

P{a, a') = , if a' = a\ieA (5.2) 

P(a,a) = l-^P(a,aO (5.3) 

ieA 

and all others are zero. 



We will use the estimate given in Theorem |2r^ to analyze this dynamics in a finite 
volume A, under a positive magnetic field, in the limit when /3 | cx). 
Let —1 and +1 be the configurations full of minuses or full of pluses, respectively. 
We will show in Lemma ^]2| that { — 1, +1} is a set metastable states. Apart from 
this characterization, we will only use the description of the energy landscape given 
in [NS], [AC,BC] and [N] in dimension 2, 3 or larger, respectively. We will show that 
the methods of Theorem p.l| allow to improve the known estimates without requiring 
further analysis of the energy landscape. In dimension 2 and 3, the improvement 
amounts to the computation of the exact (including the pre-factor C_i^+i) asymptotic 
value of the expected transition time t^ needed to reach +1 starting from —1 (that 
by Theorem p.l| is the inverse of the spectral gap of P). In higher dimension, where 
our knowledge of the energy landscape is not so detailed, we cannot compute the 
pre-factor but we show that it is a constant independent of (3, while previous results 
only gave sub-exponential bounds. 
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We remark that, unlike the exponential factor exp(— /3r(— 1)) (that only depends 
on the graph structure), the pre-factor C_i,+i is related to the particular Glauber 
dynamics we choose. 
We consider < h < 1 and assume that \fk < d — I, ^ is not an integer. 

Weset£,:=[^l. 

A (i-dimensional parallelepiped with all sides of length £ — 1 or £ is called quasi- cube 

in dimension d with maximal side-length i. 

Given a d- dimensional parallelepiped (li x I2 . . . x l^) and a (rf — l)-dimensional 

configuration r^'^^^ G ^'^'^^'^■■■^'''^ let us consider a configuration where the sites in the 

parallelepiped as well as the sites of the form /i + l, ^2, • • • , id where ri'^~^{i2, ■ ■ ■ yid) = 

+1 have plus spin and all other sites have minus spins. For such a configuration, 

as well as for all its rotations and translations, we say that r]'^^^ is attached to the 

parallelepiped. 

Following [N], we introduce a set B'^{v) in a recursive way: let B^{v) be the set of 

configurations where the pluses form a slab with volume v, B'^{v) is defined as the set 

of all configurations with volume v in the form of the rf- dimensional quasi-cube with 

maximal volume v' < v with a (rf — l)-dimensional configuration t] G B'^^^{v — v') 

attached to one of its largest faces. Heuristically, these configurations are as close as 

possible to a cube. It is easy to see that the energy is constant in every set B'^{v); 

we will denote this energy by H{B'^{y)). 

We make use of the following Theorem 

Theorem 5.1. (Theorem 3 in [N]) 

In the whole d-dimensional lattice l/, B'^{v) is a subset of the the minimizer of the 

Hamiltonian in the manifold with volume v. 

This result can be transported to the torus A(L) only for sufficiently small values 
of v/ L'^. For large values of f , the boundary conditions affect the shape of the 
minimizing configurations. Let m := min |f > l;H{B^{v)) < //(— 1)}. We take 
L so large that the configurations in B'^{v) are minima of the energy among the 
configurations with volume v for all v < m. Clearly, such an L exists, since the 
configurations winding around the torus have at least magnetization L. 
We define the set B'^ of the candidate saddles in dimension rf in a recursive way: 

1. in one dimension it is the set of configurations consisting of a single plus spin 
in the sea of minuses. 

2. in dimension d it is the set of configurations in which the pluses form a quasi- 
cube with one side of length id — ^ and all other sides of length id with a 
{d — l)-dimensional candidate saddle attached on one of the squared {d — 1)- 
dimensional faces. 

Notice that H{B'^) = max^<id H{B'^{v)). 
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Clearly, all the candidate saddles have the same volume f^ and are in B'^{v^). More- 
over, each rj e B'^{v) is connected to a configuration in B'^{v + 1) and to one config- 
uration in ;B'^(f — 1). 

Hence, the candidate saddles are saddles between —1 and +1 since from any candi- 
date saddle there exist a path leading to —1 and a path leading to +1 both reaching 
their maximum energy in the starting point. 
The following lemma was communicated to us by E. Olivieri [dHNOS]. 

Lemma 5.2. The set M. := { — 1, +1} is a metastahle set. 

Proof. We have to show that for some 5 > 0, for any a ^ —1, r(cr) < r(— 1), i.e. for 
any a ^ { — 1,+1}, there exists a configuration a' such that 

1. H{a') < H{a) - 6 

2. H{a, a') - H{a) < H{-1, +1) - H{-1) - 6. 

For 1] & Q, let \ri\ and p{ri) be the number of pluses and the number of pairs of nearest 
neighbors with different spin (namely, the perimeter, or cardinality of the contour), 
respectively. It is a well known fact that the Hamiltonian of the Ising model can be 

written as 

Hirj) = pir])-h\r]\ + Hi-l) (5.4) 

Let m := min {k > 1; 3r] with \ri\ = k and H{ri) < /J(— 1)}. Let u : —1 -^ +1 be a 
monotone one-dimensional subgraph such that ujk G B'^{k) that reaches its maximal 
energy in iS(— 1,+1). Clearly, H{ujm) < H{—1). Let a U rj denote the configu- 
ration where {a U rj) (x) := a{x) V t]{x) and a (1 f] denote the configuration where 
(a n 1]) (x) := cr(x) A ri{x). 
A direct computation shows that 

p{a) + p{r])>p{aUri) + p{anr]). (5.5) 

Since a is neither +1 nor —1, there exists at least one pair of nearest neighbour sites 
i,j such that a{i) = —1 and a{j) = +1. By translation invariance we may assume 
that the first occupied site in the sequence Uk is i and the second is j. Thus in the 
first step, aHuJi j^ —1 and H{a Uui) —H{a) < H{ui) — H{—1), while in the second 
step a U LJ2 = cr Li Ui, so that \a H u!k\ < k for all k > 2. We choose a' = aU Um- 
In order to prove point 2. we notice that for k < m, 

H{aUuJk) - H{a) = p{cr U Uk) - p{cr) - h{\(r U u^l - \cr\) 

< p{uk) - picrDUk) -h{\uJk\-\(ynuJk\) 

= H{uJk)~H{anuJk)<H{uJk)-H{-l), (5.6) 

since by definition, \a r\ujk\ < m. 

For k = m, from ( ^.6|) we get point 1. since 



H{aUuj„,)- H{a) <H{uj„,)- H{-1)<0 (5.7) 
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By putting together (ISTTI) and ( |5.6| ), we see that the energy of the configuration a' 
is lower than the energy of a and that the maximal energy in the one-dimensional 
subgraph cr U tUfc : a ^ a' is lower than if (— 1, +1). D 

In the next Theorem, we use the results of [NS] and [BC] to describe the energy 
landscape. 

Theorem 5.3. (from [NS] and [BC]). In dimension 2 and 3, 

1. The set of saddles between —1 and +1 coincides with the set of candidate sad- 
dles. 

2. //2(c/— l)//i is not an integer, the saddles are simple. 



Conjecture 5.4. Theorem \5. 5| holds in any dimension. 



We define Dj^ as the set of all states that are larger than a candidate saddle (namely 
of the form cr U r/ 7^ a for some candidate saddle a) . By the same procedure of the 
proof of Lemma ^.2| , we can easily see that all configurations in D+ have the saddle 
with +1 below the saddle with -1. We set Z := dD+ and £>_ := fi \ {D+ U Z). For 
any configuration in Z, we define pz := Tlyif^D. Pi^iV') and "pz '■= J2x'&d P{z,x'). 
Notice that in a Metropolis dynamics, all the transitions associated with an energy 
gain have probability 1/|A|. Hence, pz (resp. pz) is exponentially close to the number 
of nearest neighbor of z that are smaller (resp. larger) than z. Notice that Z is the 
outer boundary of the connected set D_ that contains —1. A direct computation 
shows that the set of candidate saddles coincides with the set of minima in Z. Hence, 
in dimensions 2 and 3, and whenever the Conjecture 5.4 holds, C_i^+i is exponentially 
close to the factor 



z€S 



^^ Pz+Pz 



1,+1 



computed with respect to Z. 
Let 



T,:=2J2 {{kil-' -{k~ l)it-') -h{il- £^l)) +2-h (5.9) 



k=2 



be the activation energy of the candidate saddle in dimension d. 

Theorem 5.5. For the Ising model on a (sufficiently large) d — dimensional torus 
A(/), in dimension d> ?>, there exists a constant Cd such that 

Er- = Qe^^'*(l + o{e'^^)) (5.10) 
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If Conjecture f)A\ holds and sk/h is not an integer for all A; = 1, . . . , rf — 1, then the 
pre-factor Cd is equal to 



^'y(i-|)n(^'^+i-^'^ + i) j (5-11) 



Lemma 5.6. The number of candidate saddles in dimension d contained in a d- 
dimensional cube of side-length I > Id is 

Nd{i) = 2''-'d\ {i-id+ 1)' n (^'^+1 - ^'^ + 1)' (5-12) 

fc=l 

All the candidate saddles have p = l^'^, while p can take the value l~'^ or 21^'^. The 
fraction of candidate saddles with p = /"'^ is j-, independently of d and I. 

Proof hetMd:=Md{Q. 

The key observation is that the pluses of a candidate saddle are contained in exactly 

one cube of side-length Cd- 

Ad — dimensional cube of side-length / > id contains {I — id + ^Y of such cubes. 

Rence,Ard{l) = {l-id + iyAfd. 

Given a cube of side- length /, there are 2d possible choices for the incomplete face 
and Md-i{id) ways to arrange the {d— l)-diniensional candidate droplet on this face. 
Hence, Md = 2dMd-i{id) = 2d{id - id-i + ly-^Ki-i- 
Since A/i = 1, a simple calculation gives ( p.l2| ). 

The computation of the number md{l) of candidate saddles with p = l~'^ is very 
similar: In dimension two, ^712(^2) = 8 i.e. the number of configurations made of a 
quasi-square plus a protuberance at one end of one of the longest sides. All other 
candidate saddles have p = 2l~'^, since there are two neighbors of the protuberance 
that can be occupied. All candidate saddles have p = l~'^, since we can void the 
occupied site and reach a quasi-cube in D_. In general, for d > 1, the only sites 
with d plus-neighbors are in an incomplete face of the d- dimensional critical cube. 
Hence, md{id) is equal the number of (c?— l)-dimensional critical squares on the faces 
of the d- dimensional critical cube times md-i{id) namely, md{id) = 2d{id — id-i + 
lf-^md-i{id-i)- On the other hand, md{l) = 2d{l -id + lYmd{id)- Thus, the ratio 
i^diJ')/Afd{l) does not depend on c? or on / and is equal to 2/i2. D 



Corollary 5.7. The number of candidate saddles in dimension d contained in a 
d-dimensional torus of side-length I > id is 

d-l 

Md{l) = 2^-'d\ /^ n (4+1 - 4 + 1)' (5.13) 

fc=i 

= d\2-^h — —{l + o{h)) (5.14) 
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All the candidate saddles have p = l^'^, while p can take the value Z""' or 2/"°'. The 
fraction of candidate saddles with p = l~'^ is f-, independently of d and I. 

Proof. The result is a straightforward consequence of lemma |5]^ and of the fact that 



the number of (i-dimensional cubes of side-length id that can be put into the torus 

is l'^. 

The estimate in ( ^.14|) comes from the approximation i^ — id-i + ^ ~ f (^ + ^(^)) 



and hence ( |5.14|) . D 



Proof of Theorems \2.2^ and \5. d^ : The results of Theorem |2]^ are straightforward con- 
sequences of Theorem p.l| , Lemma |5.3| , Lemma 5^, and Corollary 5]7. In higher 



dimension, the corresponding result comes from Theorem p.l| and Lemma |5.2|. If 



conjecture |5]^ holds, Corollary ^?^ gives the estimate in p.ll| . D 



In conclusion, let us notice that the form of the quantities P (r^ < rj^) in the case of 
the Metropolis dynamics may offer an interpretation in terms of "free energy of the 
set of saddle points" . Indeed, every point in Z gives a contribution to the pre-factor 
Cx^y that does not depend on (3 and can be bounded by a constant c. With the 



arguments in the proof of Lemma 3.2, we get 



P (r^ < T^) ~ cMe-^^ = cexp (-/5(r - T log AT)) , (5.15) 

where M is the number of saddles between x and y and T = f3~^. The logarithm 
of M can be interpreted as an entropy. This interpretation could be related to the 
results by Schonmann and Shlosman (see [ScSh]) on the connections between Wulff 
droplets and the metastable relaxation of kinetic Ising model. 

6. Appendix. 

In this appendix we briefly explain how our general approach can be generalized 
to situations than the saddles are more complicated when the isolated single points 
assumed in Section 2. The point we want to make is that in such a case it is still 
possible to localize the problem to the understanding of the neighborhood of the 
saddle points and to thus reduce the analysis of the capacities to a 'local' variational 
problem. Let us consider a situation when in the computation of a transition from 
X to y we encounter a set of saddles S^^y that can be decomposed into a collection of 
disconnected subsets iS^^\ k = 1, . . . , L. By definition, it must be true that each of 
the sets S^''^ is connected to two subsets 7?,*^^^ and M^'^'^ of D^. and Dy, respectively. 
Let us define 

C{k) := J2 e-^(^(-)-^«)p (r^(., < r>,) (6.1) 

where P is the law of the chain where all the edges exiting from the sets S^'^^ not 
leading to M^^'^ or TZ^'^^ are cut. Note that it is not difficult to see that 
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(-(fe) _ ij^f $(/^) (g 2) 

Repeating the steps of the proof of Lemma 4.2, one obtains then that 
Lemma 6.1. In the situation described above we have that 

L 

p(^;<^:) = Ec'^'^(i+«(e"'')) (6-3) 

k=l 
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